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"'Y' \ Brueckner calculations including a microscopic three-body force have been 

Q ' extended to isospin asymmetric nuclear matter. The effects of the three-body 

3 ■ force on the equation of state and on the single-particle properties of nuclear 

G ■ 

• • . matter are discussed with a view to possible applications in nuclear physics 

i and astrophysics. It is shown that, even in the presence of the three-body 

^ ' force, the empirical parabolic law of the energy per nucleon vs isospin asym- 

' metry (j = {N — Z) j A is fulfilled in the whole asymmetry range < /? < 1 

up to high densities. The three-body force provides a strong enhancement 
of symmetry energy increasing with the density in good agreement with rel- 
ativistic approaches. The Lane's assumption that proton and neutron mean 
fields linearly vary vs the isospin parameter is violated at high density in the 
presence of the three-body force. Instead the momentum dependence of the 
mean fields is rather insensitive to three body force which brings about a linear 
isospin deviation of the neutron and proton effective masses. The isospin ef- 
fects on multifragmentation events and collective flows in heavy-ion collisions 
are briefly discussed along with the conditions for direct URCA processes to 
occur in the neutron-star cooling. 
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1. Introduction 

The equation of state (EOS) of neutron-ricli matter is a source of important theoretical 
predictions on the properties of neutron stars, heavy-ion colhsions (HIC) and nuclei at 
the neutron drip line |jl|]. The interest has been mainly focussed on the properties of the 
symmetry energy, including its dependence upon the baryonic density. Studies of neutron- 
star cooling spin-polarized states 0, collective flows and isospin distillation in HIC 
have stressed that those phenomena are very sensitive to the values of the symmetry energy 
in the respective density domains. The EOS of isospin-asymmetric nuclear matter has been 
recently studied in the framework of the Brueckner-Bethe-Goldstone (BBG) theory JT],^. 
The convergence of the BBG hole-line expansion has been assessed in recent times with high 
accuracy P], even if the the saturation properties of nuclear matter are not reproduced. 
It is commonly recognized that the missing saturation is due to the model of nucleons as 
structureless particles interacting via a bare two-body force (2BF) and one has to introduce 
three-body forces (3BF). The first microscopic model of 3BF was the Fujita-Miyazawa model 
T0| where the isobar A(1232) is excited in a pion exchange interaction between two nucleons. 



The model was later extended to the A^(1440) Roper resonance [|ri|]. It has been also 
recognized that the main relativistic effect introduced by Dirac-Brueckner approach is the 
excitation of the negative energy states of the Fermi sea which can be described in terms of 
a STB |1^ as well. 

The above mentioned 3BF components have a strong saturating effect as already shown 
by a Brueckner-Hartree-Fock (BHF) calculation with the Paris potential as 2BF [0. This 



3BF has been re-adopted in Ref. |jT3| in combination with the Argonne Vis UHl as two-body 



component. In Ref. [15] it was discussed the EOS of pure neutron matter and a preliminary 
prediction fos given for the symmetry energy based on the shift between the energy of pure 
neutron matter and symmetric nuclear matter. But one may argue that 3BF could modify 
the law fulfilled by the binding energy of asymmetric nuclear matter specially at high 
density and higher order terms in the (3 expansion could significantly limit the role played 
by the symmetry energy in describing the isospin effects. All that required to to extend the 
calculation of the EOS to the full range < /? < 1 in a baryonic density range high enough. 
In this paper we present the results of such a calculation. After describing the model of TBF 
(Sec. 2) we will focus on the 3BF effects on binding energy, symmetry energy and single- 
particle properties of asymmetric nuclear matter (Sec. 3). Then a few direct applications in 
the neutron star cooling and in the isospin properties of HIC will be shortly discussed along 
with a comparison with other approaches (Sec. 4). 
2. Isospin dependent 3BF 

The microscopic 3BF adopted in the present calculation is based on the meson-exchange 
current approach.lt is described in full detail in Ref. [11]. The new meson parameters 
calculated to meet the self-consistent requirement with the adopted AVis 2BF are reported 
in Ref. [13]. 



It has been shown [|ll| , p!3| that the main contributions to 3BF rise from the two-meson 
exchange part of the NN interaction medium modified by the intermediate virtual excitation 
of nucleon resonances (isobar A(1232) and Roper A^(1440)) (Fig. la) and from the two-meson 
exchange part with nucleon-antinucleon virtual excitation (Fig. Id). This latter contains the 



relativistic effects associated to the dressed spinors of the Dirac-Brueckner approach [16 



The terms associated to the non-linear vr-nucleon coupling required by the chiral symmetry 
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T% play a minor role specially above the saturation density, where heavy mesons {a and u) 
are dominating over the 27r-3BF (Fig. lb, diagrams without heavy mesons). The contribution 
due to meson-meson coupling (Fig. Ic) is also negligible [|1T| and was not considered in the 



present calculation. In the context of non relativistic approaches phenomenological 3BF 



have also been introduced |T8| in the calculations of nuclear matter. They are modelled on 
the saturation properties of nuclear matter and/or binding energies of light nuclei. A good 
agreement between microscopic and phenomenological 3BF is found only when this latter is 



treated within the Brueckner approach 



The general Brueckner formalism for asymmetric nuclear matter (ANM) with 2BF's 
is described in Ref. [8]. The rigorous procedure to include 3BF would require to solve 
the Bethe-Faddeev equation with a 3BF much the same as already done to calculate the 
three-body clusters in the BBG hole-line expansion 0. But at this moment this appears a 
formidable task and we are content to follow a simplified procedure based on converting the 
3BF into an effective two-body force via a suitable integration over the third-nucleon degrees 
of freedom. The integral is weighted over the correlation function of the third nucleon with 
the two others [|Tl|,|T^. The effective 2BF for asymmetric nuclear matter is defined as 

<m\Vr^'\n%'> = ^EEE / df3df3'0:(r3f3')(l-r7n,r3(ri3))(l-r].,,.3(^)) 

X W-i{r(r^r'.^rxr2r-i)K{T'?,'r-^{\ - r^n.rg (?^i3))(l - ^?r2,T3(r23)) 

The function ri^^T-^ir^ is the average over spin and momenta in the Fermi sea of the defect 
function, of which only the most important partial wave components have been included, 
i.e. the ^S'o and ^5*1 partial waves. 

The transformation of the 3BF to an effective 2BF entails a selfconsistent coupling be- 
tween 3BF and Brueckner procedure of solving the Brueckner-Bethe-Goldstone equations. 
One first calculates the correlation function with only the 2BF and then builds up the effec- 
tive 3BF which in turn is added to the 2BF, and again calculates the correlation function 
and so on up to the convergence is reached. As previously mentioned the bare 2BF adopted 
in the calculations was the charge dependent AV^ig. The partial wave expansion of full 
interaction has been truncated at Imax = 6. 

3. Numerical results 

The EOS of ANM has been calculated spanning the whole asymmetry range with a 
step-size Aj3 = .2 and a density domain up to 0.45fm~'^. The case of symmetric nuclear 
matter (/3 = 0) is discussed in Ref. [15] (the saturation properties are reported in Table 1). 
The results for ANM are displayed in Fig. 2 for both cases with (left panel) and without 
(right panel) the 3BF. In the figure energy shift for asymmetric-to-symmetric nuclear matter 
is reported versus The individual runs (symbols in the Figure) are depicted along 
with their linear fits (solid lines) performed with only the first three values of asymmetry 
parameter f3. From the comparison between the two sets of calculations one may notice 
that, despite its strongly density dependent repulsive effect, the 3BF does not violate the 
law fulfilled already with only 2BF. The rather good agreement between the symbols and 
the corresponding lines (maximum deviation is 6% ) indicates the high quality of the law 
up to the largest densities. This is a quite astonishing result because the 3BF introduces a 
strong density and isospin dependence, making the nucleon-nucleon (NN) interaction quite 
different from the pure 2BF. On the other hand, it is a quite desirable result indeed for 
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two reasons. First, this (in agreement with previous studies using alone charge 



independent 2BF) provides a strong support to the apphcations based on the EOS of AMN 
extracted from the empirical /5^-law. Second, it imposes also strict theoretical constrains on 
the phenomenological nuclear forces when extended to AMN. We have in mind the Skyrme 
forces which have been fit around the saturation point of symmetric nuclear matter and 
include an effective density-dependent term to simulate the effects of 3BF. 

In Table 1 the saturation properties of ANM are shown with and without 3BF. One 
may just notice that the compression modulus results smaller with 3BF despite the strong 
enhancement of the curvature of the EOS. This is due to the fact that it is also proportional 
to the square of the saturation density, which turns out to be very much reduced. 

The symmetry energy is defined as 



Esym (p) — 2 



(2) 

/3=0 



Due to the simple /3^-law the symmetry energy can be equivalently calculated as the differ- 
ence between the energy per nucleon of pure neutron matter and symmetric nuclear matter, 
i.e., 

EsyUp) = Ea{p,1)- Ea{p, 0). (3) 

Fig. 3 shows the effect of the 3BF on the symmetry energy in the density domain considered 
in this study. At the saturation density the two values do not significantly differ: 30.71 MeV 
(3BF included) and 29.28 MeV (no 3BF). Both are in good agreement with the empirical 
value 30. ± 4 MeV extracted from the nuclear mass table ||2^. Above po the 3BF gives a 
strong enhancement of symmetry energy since it is strongly repulsive at high density. Both 
curves, with and without the 3BF, have been parametrized by simple power laws as follows, 

• BHF with pure AVis 2BF 

E,y^ = 30.7m°-^« (4) 

• BHF using AVis plus the 3BF 

Esym = 30.71m°-6 u < 1 

= 30.71 + 18.42(m-1) + 9(m- 1)2 M> 1 ^^"^ 



where u = p/po, and po = 0.17 fm~^ the empirical saturation density. The above simple 
relations, plotted in Fig. 3 (right panel), may be useful in HIC simulations. 

In Fig. 3 (left panel) Egym vs density is compared with other approaches. Despite the 



overall agreement with the predictions of both relativistic mean-field (RMF) theory |23 
and Dirac-Brueckner-Hartree-Fock approach (DBHF) the density dependence is found 
to be rather different. Both RMF and DBHF theories predict an almost linear variation of 
Egym VS. density. Instead, the present BHF calculation with the 3BF gives a slower varying 
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Esym VS. p at relatively low densities, say from p ~ O.OSfm ^ to p ~ 0.25fm ^. Such a 



■'sym 

-3 



slow variation vs. density is also found in a density region up to p ~ 0.2fm~ in Ref. ||24 
from the three-loop approximation of chiral perturbation theory. In this density region, 
the shape of Egym plays an important role as to the study of the isospin effects in HIC 
at intermediate energy as discussed later. On the contrary, in the relatively high density 
domain, i.e., p > 0.3fm~^, the present calculation predicts a steeper density dependence of 
Esym than DBHF and RMF. 

Due to the isospin effect the proton and neutron single-particle potentials are different. 
As discussed in Ref. ^ where only the charge-independent AV14 2BF was used, the attractive 
T = SD channel contribution in the two-body NN interaction mainly drives the isospin 
dependence of the proton and neutron mean fields. As a result, the proton mean field 
becomes more attractive while the neutron one more repulsive as increasing asymmetry. 
In the BHF approximation both mean fields vary linearly with the asymmetry parameter 
f3, which is in keeping with the dependence of the energy per particle. A potential 



linearly dependent on P was introduced phenomenologically long ago |25| and is referred 



to as Lane potential. Due to the importance of the single-particle properties such as mean 



field and effective mass in HIC physics p6| , |27| , it is of some interest to explore the effect of 



the 3BF on such quantities. The proton and neutron mean fields are shown as a function 
of momentum k at different asymmetries (3 = 0, 0.2, 0.4, 0.6, 0.8 in Fig. 4 for the saturation 
density p = O.lTfm^'^ and in Fig. 5 for high density p = 0.34fm~^. In both figures the left 
panels display the results using AVis plus the 3BF and the right panels those using pure 
AVis- As expected, the 3BF adds a repulsive contribution to both proton and neutron mean 
fields at all asymmetries. At relatively low density (Fig. 4), the proton mean field with 
3BF becomes more attractive, which is in agreement with the 2BF prediction. However 
the validity of the linear Lane assumption is broken by 3BF as more clearly seen in the 
left panel of Fig. 6, where the isospin variation of the proton and neutron mean fields are 
plotted at momentum k = 0. At relatively high density (Figs. 5 and 6), the 3BF force effect 
becomes much more pronounced and in fact it brings a strong deviation from the linear 
Lane assumption. The neutron mean field rises up more rapidly as compared to the results 
with pure AVis- The same happens to the proton potential which at a certain asymmetry 
becomes even more repulsive. . This remarkable result can be explained by the competition 
between the isospin dependence of the 3BF and the contribution from the attractive T = 
SD channel. As increasing isospin asymmetry the 3BF repulsion starts to compete with the 
T = SD channel 2BF attraction, and becomes the dominant one at density high enough. 

The effective mass, which is related to the non-locality (momentum dependence) of the 
neutron and proton mean fields, is defined as 



m;{k) _ k_ / dE^{k) \ ^ 
m m \ dk 



(6) 



The momentum dependence of m* is featured by a wide bump inside the Fermi sphere 
due to the high probability amplitude for particle-hole excitations near the Fermi surface 
28| . The isospin dependence of the proton and neutron effective masses at their respective 



Fermi momentum kp and kp are given in the right panel of Fig. (6). The 3BF does not 
affect the linear scissor-shaped behavior observed in the previous calculations using a pure 
2BF 0,1]. One should notice that the isospin effect on neutron and proton effective masses 
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in the Brueckner approach goes the other way around than in the RMF approach p3[ . This 
discrepancy is to be understood taking also into account the different definitions of effective 
mass in the two approaches. 

4. Summary and discussion 

In the present work the microscopic 3BF based on meson-exchange current approach has 
been extended and applied to isospin ANM in the framework of the Brueckner theory. The 
3BF effects on the isospin dependence of both the nuclear EOS and single-particle properties 
have been investigated. 

The results confirm the validity of the law for the energy per nucleon in the entire 
range of isospin asymmetry and up to high density in spite of the strong isospin and density 
dependence of the 3BF. As a consequence the isospin effects in ANM are just driven by the 
symmetry energy as in the case with only the 2BF. The vanishing of higher powers in the 

expansion also supports the simple recipe often adopted to extract the symmetry energy 
from the two limiting cases of symmetric nuclear matter and pure neutron matter . It also 
constraints theoretically the phenomenological interactions, such as the Skyrme forces which 
take into account the effect of 3BF by a density dependent term. 

As expected, the 3BF improves the saturation properties of symmetric nuclear matter by 
shifting the equilibrium density close to the empirical value. At relatively low density, the 
3BF effect on the nuclear symmetry energy is quite small. On the contrary, at high density, it 
brings a strong enhancement and consequently the symmetry energy rises with density more 
steeply than the corresponding 2BF prediction. The non linear increase of symmetry energy 



at high density has also been observed in a recent relativistic Hartree-Fock calculation ||23 



as the 'Fock' exchange effect of the non linear scalar self-interactions. But in the same 



Ref. it is also shown that at sub-nuclear densities, the Fock contribution could result 
in a softening of the symmetry potential term. 

The density dependence of the symmetry energy has been parametrized for the sake of 
application in HIC with very neutron-rich ions. It has been shown already that isospin 
fractionation Pp9| in multifragmentation events is very sensitive to the density dependence 
of Egym in the low density region. In particular, in Ref. ||3^ is shown that using an isospin 
stiff nuclear EOS with symmetry energy curvature equal to -69 MeV (very close to the 
present value of -66 MeV) leads to a remarkable value for the liquid-to-gas isospin asymmetry 
ratio which is consistent with the experimental prediction 0. An isospin scaling has been 
proposed in multifragmentation events of HIC, which turns out to be also very sensitive 
to the density dependence of Egym- Using the expanding evaporating source model and 
adopting for Egy^ the simple parametrization C ■ {p/po)"' the fragment data can be fairly 
well reproduced by taking 7 = 0.6 | JI| , which coincides with the fit of our microscopic value 
(cfr. Eq. (5)). 

The preequilibrium particle emission p2| , ^ and collective flows can instead probe 
the symmetry energy in the range of high density, say up to two times the saturation density, 
but only with an isospin stiff EOS which is in agreement with our prediction. 

In the high density range the symmetry energy is also relevant to the study of the 
neutron-star properties such as the cooling mechanism. In fact, a steep increase of Egym 
with density favours the direct URCA processes [0. In Fig. 7 the proton fraction is reported 
for /3-equilibrium nuclear matter in different approximations. The threshold of direct URCA 
processes is only crossed by the results with 3BF at a reasonable value of the neutron-star 
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density. 

As to the single-particle properties, the effect of the 3BF is to add an extra repulsion to 
both proton and neutron mean fields so that the linear Lane assumption breaks down slightly 
around the saturation density but quite strongly at high density. However, the scissor- 
shaped behavior of the proton and neutron effective masses vs. (5 remains unchanged being 
the momentum dependence of the mean fields rather insensitive to 3BF. In a calculation not 
reported here it was found that the 3BF affects only slightly the rearrangement contribution 
to the nuclear mean field and can not improve the fulfillment of Hughenoltz-Van Hove 
theorem (see Ref. [8]). This requires to go beyond BHF approximation in the expansion of 
mass operator that is a work still in progress. 
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Tab. 1 - Equilibrium density of asymmetric nuclear matter, incompressibility and energy per nucleon at 
equilibrium density corresponding to four different values of asymmetry (5. 
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FIGURES 



FIG. 1. Diagrams of the microscopic 3BF adopted in the present calculation (see Ref. [11]). 
Diagram (c) was not included. 

FIG. 2. Energy per nucleon of asymmetric nuclear matter in the range < < 1 at four 
densities as compared to the parabolic fits ( straight lines ) obtained from the first three values of 
(3 (0.0, 0.2, 0.4). Left panel : BHF predictions using AVis plus the 3BF. Right panel: BHF results 
with only pure AVi% 2BF. 

FIG. 3. Symmetry energy vs density. The left panel shows the comparison among different 
approaches: the BHF predictions with 3BF (upper solid curve) and without 3BF (lower solid curve) 
are obtained from the slopes of Fig. 2 (while symbols are the values approximated by Eq. (3)); the 
long-dash curve corresponds to the result of DBHF approach from Ref. [20]. The short dash is that 
of BHF calculation using the phenomenological Urbana 3BF in Ref. [19], where Esym is obtained 
from Eq. (3). The dotted curve is the prediction of relativistic Hartree-Fock approach taken from 
Ref. [23]. Right panel shows simple parametrizations of the present results with (solid line) and 
without (dashed hue) 3BF. 

FIG. 4. Proton and neutron mean fields in asymmetric nuclear matter at p = O.lT/m""' for 
five different asymmetries. Left part shows the proton mean-field (upper panel) and the neutron 
one (lower panel), respectively, vs. momentum using AV\^ plus the 3BF. Right part shows the 
corresponding results without the 3BF. 

FIG. 5. The same as in Fig. 4 for p = 0.34fm^'^. 

FIG. 6. Left panel: Isospin variation of proton and neutron mean fields at a fixed momentum 
A; = for two densities p = 0.17 and 0.34fm~^. Right panel: isospin dependence of proton and 
neutron eff'ective masses calculated at their respective Fermi momenta. The results in both panels 
are calculated using AV\s plus the 3BF. 

FIG. 7. Proton fraction in /^-equilibrium nuclear matter calculated in the BHF approximation 
with 2BF (^Vis) and 2BF plus 3BF (^Vis + 3BF) in comparison with a variational calculation 
using Urbana 2BF and 3BF {UVu + UVII). The dashed orizontal line is the threshold for direct 
URCA processes estimated in Ref. [34] 
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